Thermalization dynamics close to a quantum phase transition 
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We investigate the dissipative dynamics of a quantum critical system in contact with a thermal 
bath. In analogy with the standard protocol employed to analyze aging, we study the response of 
a system to a sudden change of the bath temperature. The specific example of the XY model in a 
transverse magnetic field whose spins are locally coupled to a set of bosonic baths is considered. The 
peculiar nature of the dynamics is encoded in the correlations developing out of equilibrium. By 
means of a kinetic equation we analyze the spin-spin correlations and block correlations. We identify 
some universal features in the out-of-equilibrium dynamics. Two distinct regimes, characterized by 
different time and length scales, emerge. During the initial transient the dynamics is characterized by 
the same critical exponents as those of the equilibrium quantum phase transition and resembles the 
dynamics of thermal phase transitions. At long times equilibrium is reached through the propagation 
along the chain of a thermal front in a manner similar to the classical Glauber dynamics. 



Introduction. Understanding the out-of-equilibrium 
dynamics of quantum many-body systems is a central 
issue of modern condensed matter physics from both a 
fundamental and an applicative point of view. Theoret- 
ical interest on these problems traces back to the stud- 
ies of irreversibility in non equilibrium thermodynamics. 
In quantum systems the interplay between phase coher- 
ence, strong interactions, and low dimensionality may 
result in surprising dynamical behaviors [l[. Remark- 
ably, this kind of issues can be explored experimentally at 
the quantum level by realizing highly controllable quan- 
tum many-particle systems. In this sense, cold atoms in 
optical lattices are the paradigmatic example of an in- 
teracting system where the interaction strength and the 
geometrical settings can be fine tuned. The engineered 
Hamiltonians can mimic condensed matter systems Q 
and also provide feasible tools to investigate many in- 
teresting issues in non-equilibrium statistical mechanics. 
0. Similarly, arrays of coupled microcavities have been 
shown to have the potential to act as simulators of quan- 
tum many-body dynamics, with characteristics comple- 
mentary to those of optical lattices [1] . 

In this context it is desirable to consider simple but il- 
lustrative enough situations. An important problem that 
has been studied is the response of a system to a strong 
perturbation driving it out of equilibrium. Paradigmat- 
ically, this occurs when a given system is forced out of 
equilibrium by a sudden change of a control parameter 
{quantum quench) [f|. The richness and complexity of 
this problem is ultimately related to the nonlocal char- 
acter of the correlations developing the system temporal 
evolution. An elegant picture of the phenomenon has 
been provided by Calabrese and Cardy [H, according to 
which correlations among spatially distant parts of the 
system are established because, from each point, pairs 
of entangled quasiparticles are emitted and propagate 
ballistically in opposite directions with a characteristic 



velocity. This picture, which allows to describe the dy- 
namics of correlation functions and block entropy, has 
been tested in different contexts (see |7[ and references 
therein) . 

The above results concern the idealized situation of 
the unitary dynamics of a closed system. A central prob- 
lem then, is to understand the influence of an external 
environment on the dynamics. This issue, originally ex- 
plored in the context of mean- field spin glasses [1, 0| , has 
very recently obtained some attention in the context of 
the adiabatic dynamics of open quantum critical systems 
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In this Letter we study the problem of a thermal 
quench, i.e. the response of the system to a sudden 
change of the temperature keeping the Hamiltonian pa- 
rameters fixed. This is the purely dissipative analog of 
a quantum quench [l5|. When the temperature of the 
bath in contact with the system is changed, the sys- 
tem is forced to relax towards the new thermal equilib- 
rium state. We will be focusing on such a thermalization 
process close to a quantum phase transition. Some key 
questions we wish to address are: How is thermalization 
achieved microscopically? What are the characteristic 
time and length scales emerging from the out-of equilib- 
rium dynamics? And, ultimately: What are the dynam- 
ical signatures of the quantum phase transition? 

We address the questions above for a benchmark class 
of systems with a quantum critical point, namely, the 
quantum XY model. To model a thermal reservoir the 
system is coupled to a bosonic bath. As in the unitary 
case (conservative system), we will focus on the study of 
spin-spin correlation functions. A detailed study of the 
dynamics of the system will allow us to understand how 
not only quasi-equilibrium but also far-from-equilibrium 
dynamics reveals signatures of criticality. 

Model. We consider a chain of N spins with an 
anisotropic Ising interaction and subject to a transverse 
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magnetic field: 
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where a x ' v ' z are Pauli matrices. We fix the energy scale 
J = 1 and consider /i > 0. In the anisotropic case 
< 7 < 1, the magnetic field h induces a phase tran- 
sition at h c = 1 that separates a paramagnetic phase at 
h > 1 from a ferromagnetic ordered phase with (ct x ) 7^ 0; 
such a phase transition belongs to the Ising universality 
class with critical indexes v = z = 1. The Hamilto- 
nian can be diagonalized in momentum space, in terms 
of Jordan- Wigner fermions c&, as J2k>o ^I^fe^fe: where 
= ( C L c -kj, Hfc = — (cos k + h)r z +sinfc f y , where 
f a are Pauli matrices. The T — quantum phase transi- 
tion leaves an imprint at low temperatures, leading, close 
to the quantum critical point, to a crossover at temper- 
atures T ~ A with A = \h — h c \ the energy gap. In 
particular for T <C A the spin-spin correlation function 
is factorized into a quantum and a thermal term that can 
be described semiclassically in terms of quasiparticle ex- 
citations, while in the quantum critical region (T S> A) 
quasiparticle excitations no longer exist [23 ]. 

To model a thermal reservoir we consider a set of 
bosonic baths coupled locally to each spin [3] , such that 
the global Hamiltonian reads 



n 



Hp 
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where X } = \p{b\j + b P j) and H B ^Y.j^^b^. 
The system-bath coupling is chosen to have power law 
spectral densities 5TgAjg<S(£d — ujp) = 2au s exp(— oj/uj c ) 



[23J . The system-bath coupling we are considering, 
Eq. @, breaks the integrability of the model, inducing 
transitions between all energy levels, and thus complete 
relaxation. 

The quantum quench dynamics for the closed XY 
model was studied in [l6j |. It is customary to consider 
the physical system initially uncorrelated, e.g. by ap- 
plying a strong magnetic field h. After a quench of the 
magnetic field, correlations between parts of the system 
will start to develop because of the dynamics induced by 
the new Hamiltonian. Analogously, in the case of thermal 
quenches, we consider the system to be initially prepared 
in equilibrium with the bath at a very high temperature, 
again with no correlations because the density matrix of 
the system is proportional to the identity. After a quench 
of the temperature of the bath at t = 0, the system is 
forced out of equilibrium and eventually reaches a new 
stationary thermal state. In the following we investigate 
how such correlations develop and how thermal equilib- 
rium is eventually approached. All the results shown in 
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Figure 1: Correlation function C zz for a quench from T — 00 
to T = 0.1 at h = 1; red thick line on the right at large t is 
the thermal equilibrium C ZZ (R). During the initial transient 
(t <C a), C zz oc t 2 for all R, as marked by the fits on the left. 



the figures below refer to the Ising model (7=1) cou- 
pled to Ohmic baths (s = 1). However, the results stated 
in the text refer to the general case < 7 < 1, s>0. 
We discuss the spin-spin correlation function and later 
we consider the quantum mutual information. 

Spin-spin correlations. We consider the equal-time 
connected correlation function 

C zz (t, R) = (a?(t)^ +R (t)) - <a|(t)) (<r* +R (t)) . (3) 

In the case of thermal quenches the dynamics is purely 
dissipative. Since for weak coupling a <C 1 the dynam- 
ics of populations and coherences decouple, if the sys- 
tem starts in a mixed state no coherences will develop 
after the quench (this is consistent with the so called 
"secular approximation" [l8(). Therefore in this limit at 
each time the system is approximately in a statistical 
mixture of the Hamiltonian eigenvectors, i.e. a gaussian 
state. Hence, by exploiting this the correlation function 
can be expressed as C zz = |^ J2k>o( CkC ~ k ) s ^ n (kR)\ — 
J2k>o cos(fci?)(cj £ Cfe)^ . In order to evaluate the two 
point fermionic correlators we use the kinetic equation 
derived in lj, ll7[ within the weak coupling and Markov 
approximations. From the analysis of our results (Fig.[TJ 
two regimes can be outlined: right after the quench, for 
t <C a, correlations increase as C zz oc t 2 , while in the 
opposite limit, at times t 3> a, the system is close to 
thermal equilibrium. During the initial transient, C zz 
reaches for far distant spins {R 3> 1) values greater 
than those of thermal equilibrium, C zz (R) > C*^ (R) oc 
exp(— R/Q, with £ the thermal correlation length [22|. 
Thus, the crossover to the long-time regime displays a 
non-monotonous behavior as a function of time, so that 
that C zz increases up to a maximum value and then re- 
laxes to the thermal equilibrium value (Fig. Q]) . 

Let us analyze first the initial transient. We observe 
that, for noncritical values of the magnetic field (/i ^ 1), 



3 




Figure 2: Initial transient: snapshots of C zz (left) and d^T 
(right) at a fixed t/a <C 1 after a quench from T — oo to 
T = 0.1 and for h = 0.8, 0.9, 0.95, 0.975, 1 (from bottom 
to top). The spikes relative to C zz in the left panel mark the 
distances ft (ft) at which C zz changes sign. Dashed lines are 
plotted for comparison. Inset: ft as a function of \h — h c \; 
for 8l1, ft is calculated as the maximum of d\T which marks 
the crossover between i" 1 and L -4 scaling. 



C zz changes its sign at a certain distance ft(/i) such that 
C Z z ;s for R^ f t (see Fig. [2]). That distance marks the 
crossover between two power-law behaviors with different 
exponent, respectively R~ 4 and R~ 2 , and close to the 
critical point it diverges as 



ft oc \h - h c \~ 



(4) 



Collecting all the above results, we find that the long R 
behavior close to the critical point is described by 



Czz <X t l 



R # < ft 
R- A R > ft 



(5) 



At equilibrium, close to the phase transition, the 
correlation length f marks the crossover between the 
exponential decay for R 3> f and the critical power-law 
for R <C f . Similarly, in the present non-equilibrium 
case ^ can be interpreted as an effective crossover 
scale between two power-law regimes. Eqs. Q and 
([5|) are independent of the specific value of the final 
temperature at which the system is quenched and of 
the specific exponent s of the bath spectral function. 
Moreover, they are robust within the range < 7 < 1 
in which the system belongs to the Ising universality 
class. Remarkably, although in this regime the system is 
far from equilibrium, Eqs. (j4|) and (O are characterized 
by the equilibrium critical indexes: f oc \h — /i c | _1 and 
CzzOiR- 2 , R «f. 

We now focus on the long time regime. The analysis 
of Fig. [3] indicates that, at a given time, the correlation 
function is thermalized up to a distance Rth- This ther- 
mal front exists because the long-distance correlations 
are dominated by the slowly relaxing low-energy modes. 
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Figure 3: Thermalization of the spin-spin correlation function 
Cz Z after a quench from T = 00 to T = 0.1 at h = 1. Left: 
snapshots of C ZZ {R) at t/a — 10, 20, 30, 40 from top to bot- 
tom; thick red line is the (exponential) thermal equilibrium 
C*z- At a given time after the quench C zz is thermalized up to 
a distance Rth that increases with time. Right: corresponding 
time dependence of i? t h; the linear fit gives v t h — 0.32. 
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Figure 4: Thermal front velocity -u t h as a function of 1/T. 
From top to bottom h=l, 0.9, 1.2 (so that A = 0, 0.1, 0.2). 
Lines are the fit Ta(l + b£) e~ A/T with a = 3.3, b = 0.9 for 
the specific case 7=1, 
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The front is found to propagate ballistically with a speed 
Vth that is a function of T and h. In particular, as shown 
in Fig. [4l the velocity scales as 



T 



Vth oc 



T > A 

A/T T<<A 



(6) 



where A = \h— 1| is the energy gap of the XY model. 

Block correlations. We now analyse how correlations 
between a block of spins and the rest of the chain develop 
after a thermal quench. In order to quantify such a cor- 
relation, we use a tool originally developed in the context 
of quantum information theory. For a certain bipartition 
of the system into two blocks of L and N — L spins, the 
mutual information is defined as 



2(L) = S(p L ) + S(pn-l) - S(pn) 



(J) 



where S(p) — — Tr(plogp) and pn is the density matrix 
of the entire system. The mutual information measures 
the correlations between the two blocks of L and N — L 
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spins 20]. For the XY model T{L) is known to diverge 
logarithmically as a function of L at the critical point, 
while it saturates for noncritical values. In the following 
we concentrate on the derivative dj,T, which measures 
the sensitivity of the correlations to the block size. It 
is useful to study d^I because, at equilibrium, it shows 
features similar to the spin-spin correlation function: it 
scales as drX oc L^ 1 at the critical point, while it decays 
exponentially for noncritical values. 

Equation (|7|) can be computed in terms of two-point 
fermionic correlators (obtained b y s olving the kinetic 
equation) using the results of Ref. [2l( . We study for the 
mutual information the same setting of thermal quenches 
investigated above for the spin-spin correlation function. 
Remarkably the scenario emerging is very similar to that 
depicted in Fig. [TJ There are two regimes: an initial 
transient governed by 

OlI oc t < . 8 
\l- 4 L>& V ; 

with the same correlation length (see Fig. In the 
quasi-equilibrium regime at long times, 8lX exhibits a 
thermal front propagation similar to that shown in Fig. 
|3]and with the same velocity found for C zz (see Fig. |4|). 

Conclusions. We have analyzed the dynamics of spin- 
spin and block correlation functions following a sudden 
cooling of the bath coupled to a quantum system. For 
both quantities we find that the dynamics displays two 
regimes: at short times the correlations develop accord- 
ing to © and ((8]), while at long times a well defined 
thermal front propagates along the chain with velocity 
the latter being sensitive to the critical properties 
of the system. We remark that the system does not ex- 
hibit aging because it is quenched away from the critical 
point. Nevertheless in its early stages relaxation does 
show critical features analogously to those of thermal 
phase transitions. In particular, for systems quenched 
at the critical temperature, it is known that equal-time 
two-point correlation function scales, during the initial 
transient, as a power law both in time oc t a (here a = 2) 
and in space R- d + 2 -n ( m our case 

R- 2 ) .15]. Besides, 
we point out that the scaling of the thermalization veloc- 
ity w t h oc exp(— A/T), which we find in the semiclassical 
regions T A, holds also in the classical Ising model 
within the Glauber dynamics [3] ■ This similarity can be 
ascribed to the fact that the system-bath coupling gener- 
ates incoherent relaxation without conserving the order 
parameter, as happens in the phenomenological Glauber 
model for the dynamics of the classical spins. 
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